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ABSTRACT 



This is a continuation of our earlier work where we constructed a phe- 
nomenologically motivated effective action of the boundary gauge theory 
at finite temperature and finite gauge couphng on S"^ x . In this paper, 
we argue that this effective action quahtatively reproduces the gauge the- 
ory representing various bulk phases of R-charged black hole with Gauss- 
Bonnet correction. We analyze the system both in canonical and grand 
canonical ensemble. 



June 2007 



^e-mail: tanay@iopb.res.in 

^e-mail: niukherji@iopb.res.in 

■^e-mail: subirQiopb.res.in 

*e-mail: swarnen@iopb.res.iii 



1 Introduction and Summary 



In a recent paper [1], we initiated a study of the phase structure of gravity in anti- 
de Sitter space in presence of higher derivative corrections. Our motivation stemmed 
from the fact that due to AdS/CFT correspondence [2,3] a quantum theory of gravity 
provides natural arena for addressing issues in gauge theory and vice versa. Though 
absence of a formulation of string theory on AdS space and lack of adequate techniques 
to study strong coupling regimes of gauge theories make a quantitative comparison 
difficult, nevertheless, encouraging agreements among qualitative features have been 
gathered over the last few years. One of the remarkable steps [4] in this realm is 
the identification of the crossover from the thermal AdS phase to black hole phase, 
namely the Hawking-Pagc transition [5], on the gravity side with large deconfine- 
ment transition on the gauge theory side ^. Based on this identification, recently, a 
phenomenological matrix model was proposed [11], which belongs to the same univer- 
sahty class of M—^ supersymmetric SU (N) gauge theory on at the limit of infinite 
't Hooft coupling (A). This model, which we will call (a, 6)-model in the sequel, was 
characterized by two coefficients a and b which depend on temperature T (and A). It 
correctly reproduced the qualitative features of the phase structure of the dual theory 
on the gravity side. This model was also used to study the transition between AdS 
sohton and black hole [12] 

Following this fine of works, in [1] we considered the response of the Hawking- Page 
transition and the associated thermodynamic phase structure to higher derivative 
terms in the gravity theory, which appears as a' corrections of the underlying string 
theory. While a study with general class of higher derivative terms would be desirable, 
in [1], we restricted ourselves to the Gauss-Bonnet (GB) terms only. The reason for 
this is that the GB terms can capture non-trivial first order effects of a' correction 
while admitting explicit black hole solutions [13-19, 29] ^. It turned out that the 

'"'In many papers, it was argued that this large N transition might even occur at zero couphng [6,7] 
and it turned out that for zero couphng, the transition appeared exactly at the Hagedron temperature 
of the low temperature thermal AdS phase. As the coupling is increased the Hagedorn and the 
deconfincmcnt transitions separate out. The Hagedorn transition occurs at a higher temperature 
than that of the deconfinement transition. However it is also believed that this Hagedorn transition 
at weak coupling manifests as the deconfinement or Hawking-Page transition at strong coupling. 
For a related work see [8] . Non-perturbative eflfects near Hagedorn transition was studied in [9] (see 
also [10] 

^GB term is expected to occur only in Heterotic or K3 compactification of type HA theory but 

not in type H theories with maximally supersymmetric compactification. The lowest correction in 
type HE theory on AdS^ is of order a The thermodynamic phases of the perturbative supergravity 
solutions as well as their boundary duals have been studied by various authors [21-24]. However the 
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phase structure in GB theory depends crucially on the GB coupling. For sufficiently 
strong coupling there is a single big black hole phase which serves as a local minimum 
below a critical temperature. As the coupling reduces below a critical value, two 
additional black holes of intermediate and small sizes appear, of which, the former 
one has negative specific heat while the small hole has positive specific heat. The 
latter one is stable and corresponds to a local minimum up to a critical temperature 
beyond which it ceases to exist. Prom the associated phase diagram, we find that the 
phase structure is similar to that of a Van der Wall's gas. Along with the standard 
HP transition between the big black hole and thermal AdS, we identified one more 
phase transition where crossover of the energies of small black hole and big black hole 
occurs. 

A similar study of the various phases in the dual theory is also carried out in [1]. 
Usually, to describe gauge theory on x at zero coupling one writes down an 
effective action with vev of Polyakov loop on as the only light degree of freedom. 
Since it is hard to compute such an action at finite YM coupling, one relies on 
phenomenological model [11], that belongs to the same universality class of the gauge 
theory as described in [11] . Appealing to the universal nature of this model near 
the critical temperature, we analyzed in [1] the j dependence of the coefficients (a, b) 
by mapping the GB corrections to the gauge theory. We also suggest [1] a modified 
matrix model which, with appropriate restriction on the coefficients of the model, 
qualitatively captures the bulk properties with a' corrections. This is an extension of 
(a, 6)-model and includes higher power of density of eigenvalues of the vev of Polyakov 
loop indicating, perhaps, that other operators in matrix model become relevant. 

In this paper, we extend our study of the effects of higher curvature corrections 
to charged sector in the AdS/CPT set up with an aim to identify some universal 
features of the boundary matrix models at strong coupling. The charged sector of 
the GB action contains a Maxwell term besides the GB corrections. Maxwell term 
typically comes in type IIB theory on AdS^ x from angular momentum along 
direction, or in other words, from the 5*0(6) gauge symmetry arising from the group 
of isometrics of S^. We focus our attention to those black hole configurations which 
have equal U{1) charges for all the three commuting U{1) subgroups of S0{6). These 
black holes and their phase structures were considered in [25,26] (see also [27]). We 
study the changes of phase structures due to GB corrections in canonical and grand 
canonical ensembles. On the gauge theory side ^, in order to describe the charged 

qualitative phase structures in this case is quite similar to that of gravity with GB (a') correction. 
^We would like to emphasise that the inclusion of GB correction takes us away from type IIB 
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sector, we use the same model as in [1] except we allow the coefficients of the model 
to depend on appropriate parameters of the ensemble as well along with the temper- 
ature (T) and the t'Hooft couphng (A). In the following we summarize our results 
along with some already known facts about phase diagram in absence of GB coupling. 

Grand canonical Ensemble (Fixed Potential): 

In the grand canonical ensemble, the black hole is allowed to emit and absorb charged 
particles keeping the potential fixed till the thermal equilibrium is reached, which, 
in this case, is governed by a fixed chemical potential. Here the phase diagram is 
characterized by the chemical potential $ and the GB coupling which we call a in 
the paper. 

($ 7^ 0, a = 0) : On the gravity side, the phase diagrams have been analyzed 
in [25,26]. If the potential $ is below a critical value, various phases are similar 
to that of AdS-Schwarzschild black hole while for $ large enough, the black hole 
free energy becomes negative compared to that of the thermal AdS at any fixed 
temperature. On the gauge theory side, at zero and small A, the phase structure was 
analysed in [28,29], while for large A, a phenomenologically motivated matrix model 
can be constructed and we will have occasion to elaborate on it at a later stage of 
this paper. 

($ 7^ 0,0; 7^ 0) : This case is studied in section (2.1) where we find the critical 
value of $ depends on a. For small there are three different black hole phases; 
one of them being unstable. Identifying the rest two as a small and a big black hole, 
we find that there is a first order phase transition from the small to the big black 
hole. However, once thermal AdS is included in the phase diagram, we find both the 
small and big black hole phases are metastable at low temperature and big black hole 
becomes stable only at high temperature. In order to clearly illustrate the various 
phases in this range of $ we construct a Landau function with black hole horizon 
radius as the order parameter. When $ is above the critical value, phase diagram 
shows a single black hole phase which is stable beyond certain temperature while a 
crossover from black hole phase to thermal AdS occurs for temperature lower than 
that. If we increase $ even further, the black hole phase is always found to be stable. 
All these phases can be summarized in a (<I>^ — a) diagram; sec Figure 1. We also note 
that, in all the above cases, wherever there is Hawking-Pagc transition from AdS to 
the black hole phase, the transition temperature is found to decrease with a. 

We study the grand canonical ensemble of the dual theory in section (3.3) and 

framework. However we will assume that a version of gauge/gravity will continue to hold. 
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(3.4). Here the parameters of the matrix model depend on the chemical potential 
(/x). We find for the chemical potential less than the critical value the analysis is 
similar to that of [1] . As in [1] , the matrix model has an extra saddle point that has 
no analogue in supergravity. We interpret this saddle point with some phase in string 
theory^ Beyond the critical potential we encounter different possible situations de- 
pending on the position (expectation value of the Polyakov loop) of the extra saddle 
point. 

Canonical ensemble (Fixed Charge): 

In the canonical ensemble, the black hole is allowed to emit and absorb radiation, 
keeping the charge fixed till the thermal equilibrium is reached and the phase diagrams 
are characterized by the charge of the black hole, q and the GB coupling a. 

(g 7^ 0, a = 0) : The phase structure is discussed in [25,26] in great detail. There 
exists a critical charge above which, at all temperature, only one black hole phase 
exists. Below gc, there can be at most three black hole phases. We call them small, 
intermediate and large. While the intermediate one is unstable, the small and big 
black holes are stable. It was also noted that thermal AdS is not an admissible phase. 
When we increase the temperature, there is a crossover from a small black hole to a 
large black hole phase via a first order phase transition. 

(g 7^ 0, a 7^ 0) : This part of the analysis is given in section 2.2, where, we find 
the phase structure depends on two parameters g and a. In particular, in (g^ — a) 
plane, we identify two distinct regions (see Figure 5) where region I consists of three 
black hole phases, while in region II, only one black hole phase exists. Thermal AdS 
continues to be non-admissible phase. As before, there is a transition from small to 
big black hole at a critical temperature. This temperature decreases as we increase 
a. 

Study of the canonical ensemble of the dual theory is discussed in section (3.5). 
Since the explicit dependence of the coefficients of the matrix model on the chemical 
potential is very hard to determine at strong coupling, we assume the the zero coupling 
result is valid there or at least the universality class of the theory does not change 
once we tune up the gauge coupling. We find that this dependence is consistent with 
one of the possible scenarios. For this scenario, we write down the matrix model for 

^This could be related to some state in string theory. At this point it may be mentioned that 
appearance of string states in boundary theory also occurred in [11,30] which corresponds to the 
Gross- Witten transition [31,32] and which was identified with a crossover from supergravity black 
hole solution to string state [33] . 
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the fixed charge and find that this model correctly reproduces the phases of the black 
holes with fixed charge. 

The paper is structured as follows. We begin with the thermodynamics of charged 
sector in presence of GB couphng in section 2. Subsections (2.1) and (2.2) are devoted 
to discussion of canonical and grand canonical ensembles. The dual theory is discussed 
in section 3. We begin by reviewing zero charge sector of the dual theory. A short 
account of the zero and weak 't Hooft coupling results is given in subsection (3.1) and 
in subsection (3.2) we briefiy mention the comparison of phase structures between 
gravity and matrix model. The rest of the section 3 is devoted to our study of the 
charged sector. The grand canonical ensemble of the dual theory is considered in 
subsection (3.3) and (3.4) while the canonical ensemble is discussed in subsection 
(3.5). 

2 Gauss-Bonnet black hole with electric charge 



7 = / d"+^x^/=^[-^ - 2A + aiR" - ARabR"' + RabcdR''"^) - ^1 ■ (1) 



We start with n + 1-dimensional (n > 4) action 

R_ 

Here, a is the GB coupling. As the higher derivative corrections are expected to 
appear from the corrections in underlying string theory, we will often refer the GB 
term as a' correction in this paper. This action possesses black hole solutions which 
we call charged GB black holes [27]. These solutions have the form 



ds'^-V{r)df + :^^+r'dnl_„ (2) 



where V{r) is given by 



^2 ^2 



4q; Aam iaq 
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In the above, a — (n — 2)(n — 3)aKn-\.i, while I is related to the cosmological constant 
as f = —n{n — l)/(2«;„+iA). The parameter m is related to the ADM mass of the 
hole, M as 

M = (4) 
where Un-i is the volume of the unit (n — l)-sphere. The parameter q gives the charge 



^ 2^2(n-l)(n-2)^._ig 

W — ) [^) 
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of the electric gauge potential 



n 



1 



2(n - 2) r"-2 



(6) 



where $ is a constant which we will fix below. Denoting r+ as the largest real positive 
root of V{r), we find that the metric (3) describes a black hole with non-singular 
horizon if 

(^)r2-2 + zV">g¥. (7) 



Finally, we shall choose the gauge potential to vanish at the horizon. This fixes 
$ to be 



n — 1 



(8) 



/ 2(n-2)r+'^-2' 

This quantity is the electrostatic potential between the horizon and infinity. Asymp- 
totically, the metric (3) goes to AdS space, as in this limit. 



V{r) = 1 + 



1 



-2a 2a 

Hence we notice that the metric is real if, 

a < —. 



(1 



(9) 



(10) 



We shall restrict ourselves to a which satisfy the above bound. In this paper, we 
will primarily consider black holes in five dimensions (n = 4). However, it is easy to 
extend the results of this section to higher dimensions. 

The thermodynamic properties of the black hole will depend on whether we con- 
sider the canonical ensemble (fixed charge Q) or grand canonical ensemble (fixed 
potential $) . The equilibrium temperature T can identified from the period /3 of the 
Euclidean time of the metric (3), which in five dimensions is given by 



2nr+{rl + 2a) 



(11) 



rl + 2rl/P -qyrl' 
As it will be useful for us to write thermodynamic quantities in terms of dimen 
sionless quantities, we define 



a 



Q - 



m 



(12) 



In terms of these quantities, (11) can be expressed as 
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Figure 1: The curves in the (q;-$^) plane separating various regions with different 
behaviours of black holes. 



27r/f(f2 + 2a) 



(13) 



2.1 Grand canonical ensemble 



In the grand canonical ensemble, with fixed potential $, the free energy can be 
computed from the Euclidean continuation of the action (1). We obtain the action 
(subtracting the AdS background) as 



+ (18a - 1 + 4$V3)f^ + 3a(l - 8$V3)r^ - 



(14) 



where (5 is inverse of T expressed in terms of potential 

27r/(f2 + 2a) 



f(l -4<I)V3 + 2f2)' 



(15) 



It will be important for us to find out the number of turning points of f3{f) as we 
vary a and $. First of all, the nature of f3{f) depends crucially on the value of 
For <l>^ > 3/4, f3{f) blows up at = (4$^/3 — l)/2. Consequently, the temperature 
is zero. Following [25,26], we would like to identify this with an extremal hole. For 
f less than this value, f3{f) becomes negative. It can easily be checked that as long 
as $^ > 3/4, regardless of the value of a, there is no turning point of /3(f). If, on the 
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other hand, <l>^ = 3/4, /9 diverges at r = and goes to zero for large r. Now, to have 
turning points, dl3/df = 0. This gives 



^^'2 12 



^ (3 - 36a - 4$^ ± ^(-3 + 12a + 4$2)(_3 + I08a + 4$2)) . (I6) 



Prom here, it follows that in order to have real roots, a should not lie within the 
window 

1 / 4$S 1 / 4$S , , 

36^ 3^~~4V 3^ ^ ' 

However, it is easy to check that for o; > |(1 — '^1^2 negative while a < 
^(1 — ^), positive. Hence, we have the following picture. For $ < \/3/2, /? 

has two turning points only if 

a<—(l . (18) 

- 36^ 3 ^ ^ ^ 

The above features of j3{f) can be nicely summarized is a — a) diagram. This 

is shown in Figure 1. The region satisfying (18) is the region 1 in the figure. So, here 

f3{f) has two turning points. However, note that for $^ < 4/3 and a = 0, j3{f) has 

only one turning point at non-zero f. In the rest of the regions namely 11, 111 and IV, 
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there are no turning points of (3{f). However, as in I, in region II, f3{f) diverges at 
f = while in regions III and IV, f3{f) blows up at finite non-zero values of f. There 
are other differences in these four regions (particularly when the free energies of the 
black holes are considered). This is what we discuss in the next paragraph. Various 
representative plots of the (3 versus f for all these regions are shown in Figure 2. 

Free Energy: The free energy can be obtained from (14) as 1^ = Igc/ P- For different 
values of the parameters, W has been plotted as a function of temperature in Figure 
3. In this figure, (a) corresponds to the parameter values where we have only one 
stable black hole solution. This is also the situation in the case of (d). However, there 
is a distinct difference in their phase structures as is evident from the plots. While the 
black hole phase has lower free energy than thermal AdS in (a) for all temperatures, 
there is a Hawking- Page transition at Tc in (d). This difference can easily be located 
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in — a) diagram in Figure 1. In the region IV of this figure, black hole at T = 
has less energy than the thermal AdS and hence is stable. However, in regions II and 
III, a hole at T = is in a metastable phase while AdS is the stable one. Hence this 
hole would decay to AdS by radiating away its energy. The fine separating HI and IV 
represents hole with zero free energy at T = 0. The equation of this fine as a function 
of a and is obtained by setting W{f, a, = 0, where = 1/2(1 - 4$^/3). Note 
that this also means that on this line the Hawking-Page temperature vanishes. 

Now returning back to Figure 3(b), we see that at low temperatures there is no 
black hole phase. Two black hole phases appear as we increase the temperature. 
The small one turns out to be unstable and the large one undergoes a Hawking 
Page transition at Tc. Note that Figure 3(c) is similar to the one we found in our 
previous paper [1] (i.e for <l>^ = O.a 7^ 0). As in [1], we have therefore the following 
scenario. At low temperature, free energy has only one branch (branch I). However, 
as temperature is increased, two new branches (branch II and III) appear. Branch II 
meets branch I at a certain temperature (Ti) and they both disappear. On the other 
hand, branch III continues to decrease, cuts branch I at a particular temperature 
(T2) and becomes negative at temperature T^. These three branches represent small, 
intermediate and large black hole. Out of these three, the intermediate black hole is 
unstable with negative specific heat, while the rest are classically stable. As in [1], we 
get two first order phase transitions (HPl, HP2). HPl is a transition between small 
and large black hole at temperature T2 and the other (HP2) is the usual transition 
between AdS to large black hole as Tc. 

It is easy to construct a Landau function which represents the behaviour of the 
free energy around the critical points. Identifying f as the order parameter, this 
function can be written as 

>V(f, T) = ^ 3f^ - 47r/Tf^ + (3 - 4:^^)f^ - 2A'KalTf + 3a . (19) 

Notice that this expression reduces to the one in [34] when we set a to zero and to 
the one in [1] as we set $ to zero. It can be checked that at the saddle point of 
this function, we get the temperature (given by the inverse of the expression in (15)). 
Substituting back the temperature in W, we get the free energy W. We have plotted 
the Landau function in Figure 4. Consider (c) in Figure 4 in particular. This is for 
= 0.2,0: = 0.01. Clearly, for T < T2, the global minimum occurs for small f, 
representing a stable small black hole phase. At T = T2, small and big black hole 
co-exist. At even higher temperature, big black hole represents the stable phase. 
However, all these phases are meta-stable below T < Tc 'ii we include thermal AdS 
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Figure 4: Landau function W vs. r for different temperatures, (a), (b), (c), (d) 
correspond to values of and a of Figure(2) 

(representing the horizontal line with W = 0). Big holes then are only stable beyond 
Tc. Note, that for (a) in Figure 4, black hole is always the stable phase while for 
(6), there is a crossover from thermal AdS to black hole at Tc. Finally, Figure (d) 
is similar to (a) except that the f — hole has more energy than the thermal AdS. 
To this end, we note that since W — E — TS — ^Q, we get the energy, entropy and 
charge as 



where expressions for Q, M in terms of g, m respectively were defined earlier. It can 




(20) 
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Figure 5: The curve va. a — (f plane separating regions with various number of solu- 
tions. 

be checked that the first law of thermodynamics dE — TdS + ^dQ is satisfied. 
2.2 Canonical ensemble 

Wc now consider the system in canonical ensemble where the charge q is kept fixed. 
We first note from the expression of the temperature (11) that T is non-negative if 



When the equality is saturated, the temperature is zero and we call this an extremal 
black hole. Denoting mass and the scaled horizon radius as fhe and fg respectively, 
we see that the following relation is satisfied: 



Like in the previous case of fixed potential, we now identify the relevant regions 
in the {a-(f ) plane. The curve separating the regions for various number of positive 
solutions for f is given by the following parametric equations in f: 



> 0. 



(21) 



m, 




(22) 




(23) 



a 



5 /f2 - 3f^ 
3 V18f2 + 2 



) 
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The curve is shown in Figure 5. It can easily be checked that for any point in region 
II, there is one real positive root for f at any temperature. Below this, that is in 
region I, there is a maximum of three. Furthermore, unlike the fixed potential case, 
the vertical ^^-axis i.e. for a = we also have a maximum of three real positive 
solutions as long as < ql{= 1/135). We also notice that thermal AdS exists only 
when — 0. The corresponding (5-f plots for these regions are shown in Figure 6. 

Free Energy: We now compute the action (1) in the fixed charge ensemble. After 
properly adding a boundary charge and subtracting the contribution to the extremal 
background, we get 



_ U^fp r 2 _4 m{f-f^-2f^) 3_2 %f 

- ^ + 7^ + f2(f2 + 2a) 2''^ ~ 2f2 



K5 



(24) 



where (3 is 

^ 2nl{f' + 2a)f 

The free energy can therefore be obtained as F = Id (3. Behaviours of free energy 
for different values of (g^, a) are shown in Figure 7. We first of all note that, in 
fixed charge ensemble, black holes with negative free energies are always the stable 
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compared to thermal AdS. Secondly, in (a) and (d), we see that given any temperature, 
there is a single black hole phase, while in (6) and (c), there can atmost be three 
phases. We call them small, big and intermediate black hole phases. We find that at 
a certain temperature, which we call T later, there is a first order phase transition 
from small to big black hole phase. On the other hand, the intermediate black hole 
is an unstable phase with negative specific heat. It can be shown by comparing (b) 
and (c), that T decreases as a increases. 

Finally, the Landau function can be constructed as before. It is given by 

^(f, T) = ^ \3f^ - ATTlTf'' + 3f2 - 247r/aTf + ^ - S - ¥l ■ (26) 



K,5 



r 



It can be checked that, at the saddle point, it reproduces correct temperature T and 
the free energy F. A plot of this function for different temperatures is shown in 
Figure 8. As can be seen in (a), for high g, there is a single black hole phase for 
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Figure 8: The Landau Function JF as a function of f for different temperatures. The 
critical temperature at which there is a transition between the small black hole to 
the large black hole is shown by T in the figures. Finally, (a), (b), (c), (d) represent 
plots for values of and a of Figure 6. 

all temperatures. As we reduce q beyond certain value, two new black hole phases 
appear in (6) for a certain range of temperature. Above and below this range there is 
only one black hole solution. When the temperature is within this range, for T < T, 
the small black hole is favoured. Otherwise, big black hole is the stable one. There 
are degenerate minima at T = T representing phase co-existence. 

Now as we turn on o!, wc get (c) for low values of q, a. This is similar to (6) except 
that the critical temperature T reduces with a. Again, for large a, we get (d) which 
is qualitatively similar to that of (a) . 

Finally, since F — E — TS, we have 
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It can be checked that they satisfy the first law of thermodynamics dE = TdS + ($ — 

3 Dual Matrix Model 

Classical Type IIB supergravity is dual to jV = 4 SU (iV) gauge theory at strong 
coupling. The phases of the gauge theory in this regime of the coupling is very hard 
to analyze. Assuming the AdS/CFT correspondence to be true one can on the other 
hand get some understanding from the supergravity computations, that are more 
tractable. Following these hues, Witten has argued in [4] that the Hawking-Page 
transition in gravity corresponds to the confinement/deconfinement transition on the 
boundary. This is the philosophy that we will follow. 

In this section, using the AdS/CFT correspondence we write down the dual matrix 
model corresponding to the i?-charged black holes that were studied in the earlier 
sections. The matrix model is constructed phenomenologically so that it captures the 
phases of the black holes. The results from the weak coupling computations available 
in the literature will be used as a guideline in this prescription as was also followed 
in [1]. Wc begin with a brief review of some of the results that we will need for the 
finite temperature gauge theory at weak coupling. 

3.1 Weak coupling 

jV = 4, SU (N) gauge theory at zero and weak couplings has been analyzed by various 
authors (see for example [6,7]). For large N, when the t'Hooft coupling A = Qym^ is 
small or zero, some of the results are explicitly known. Specifically, when A = 0, it was 
shown that the boundary gauge theory at finite temperature on x undergoes 
a phase transition that can be identified as the "deconfinemcnt" transition. S"^ x 
is a compact manifold and thus allows only colour singlet states by the Gauss law 
constraint. Though non singlet states are never possible, there are various indications 
that this transition mimics the deconfinemcnt transition in gauge theories. One of 
the indications is that there is a jump of the free energy from order to order N'^ 
while another is a discrete change in expectation value of the Polyakov loop. 

Dimensional reduction of the A/" = 4 theory on x 5*^ leaves only one massless 
mode, namely, the zero mode of Aq. One can thus write down an effective action 
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by integrating out all massive modes. The resulting model with the gauge fixing 
conditions, diAi — and dfa = 0, is a zero dimensional matrix model given by, 

Z ^ J DUe^'"^^\ (28) 

where cos denotes the volume of S^. Apart from the free energy, whose N dependence 
reflects the phase of the gauge theory, one can also define a Polyakov loop that acts 
as an order parameter for the dcconfincmcnt transition. The Polyakov loop defined 
along the time circle of x is (l/A^)TrPexp(i Jq dtAo). which is nothing but 
(l/iV)Tr(f/) as can be seen from (28). It can be shown that the expectation value of 
the Polyakov loop does indeed vanish at low temperatures and picks up a non-zero 
value above some finite temperature, which may be called the Hagedorn or the 
" deconfinement" transition temperature. 

The explicit form of the partition function is given by. 



/CO -j 
dUexp -^(^") (Tr(C/")Tr(C/-")) 

^(X") = Zy(x")+^5(2;") + (-l)"+^-ZF(2;") ; X ^ 6 



(29) 



Introducing the density of eigenvalues for U and defining p = {l/N)Ti{U) one can 
write down the effective action as Seff (see [7] for details). In the large N approxi- 
mation, the various phases are given by the solutions of the saddle point equations of 
motion. In terms of the density of eigenvalues, the above transition is reflected by a 
jump of p from zero to a nonzero value below and above Th respectively, where Th 
is given by z{x) — 1. 

These features are some what modifled when a small value of the coupling A is 
turned on. It is possible to write Seff only in terms of powers of Tr(C/) by using the 
saddle point equations. An effective action containing only the quartic interactions 
may be written as ^, 



^This model is obtained by keeping tains upto ©(A^) in the effective aetion. In the large N limit 
such terms come from three loop computations. It also determines the sign of b. In [7] the phases 
for both the signs of b have been studied. An explicit computation for pure Yang Mills theory on 
a three sphere at finite temperature shows that b is positive, implying that the transition is of first 
order at weak coupling [35]. 
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Z{X,T)^ I 



dUe^p a(A,T)IV(C/)IV(C/)t + A(A,T)(lV(C/)IV(C/)t)' . (30) 



The equations of motion resulting from (30) are, 



ap + 2bp^ 



P 



1 1 




(31) 



4(1 - p) 2 



-<p<l. 



The matrix model (30) undergoes two different phase transitions as a function of 
temperature. One is a first order transition like the zero coupling case, when 6 > 0. 
The other is a third order transition for which the eigenvalue distribution goes from 



3.2 Strong coupling and comparison to gravity 

If one assumes validity of (30) in the strong coupling regime (where a and b are some 
complicated functions of A) one may hope to map these phases to those of gravity 
obtained in the supergravity approximation. It was shown that this simplified model 
indeed possesses the thermodynamic behaviour expected from gravity when a < 1 
and 6 > [11]. A new feature that is not visible in the supergravity approximation is 
the transition corresponding to the third order transition that was mentioned in the 
earlier paragraph. It was conjectured that this should correspond to black hole string 
transition [33] . It should however be noted that the matching is only qualitative and 
the vahdity of the effective action (30) in the strong coupling regime is only limited 
to the regions around the critical points. 

The phase structure gets modified once we include (a) corrections. In addition to 
the large stable black hole, there is also a small stable black hole solution along with 
an intermediate unstable one. There are two phase transitions: (i) Prom the small to 
the large black hole (HPl) at T2 (ii) Between thermal AdS and the large black hole 
(HP2) at Te. 

In [1] we studied the matrix model that captures these phases in gravity. Since 
the size of the small black hole is less than a/cZ while we include corrections only 
up to first order of a' we discussed both the following possibilities. If we do not 
trust this solution, the situation is same as that of the case without higher derivative 
corrections and the (a, b) model (30) itself serves as the dual of gravity. With the a 



the gapless phase for < p < 



to a phase with a gap for ^ < p < 1. 
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corrections in the bulk it is also possible to see the variations of a and 6 as a function 
of l/y/X. An interesting point is to note that b decreases as we decrease A. 

On the other hand, if we trust this solution the minimal action for the matrix 
model gets modified so that it can reproduce all the features of gravity near the 
critical points, including the small black hole. This modified action is given by. 



(32) 



The Aj's are functions of A and T. The (a, h) model is recovered when A^, and 
vanishes. The addition of a corrections in gravity requires higher powers of p to 
be added. In order that the model (32) reproduces all the features of gravity one 
also needs to restrict the values of AiS within certain regions. There is however one 
saddle point solution that does not have a counterpart in supergravity. This solution 
is locally unstable and can possibly serve as some stringy decay mode for the small 
black hole. In the following parts of this section we will use (32) so that it incorporates 
the phases of charged black holes as discussed in section (2). 

3.3 Matrix model with chemical potential 

Once we turn on a non-zero chemical potential the coefficients of the matrix models 
described above will depend on the chemical potential. At zero coupling, ie at A = 0, 
this dependence is easy to determine. The partition function in this case is given by 
(29) with z{x^) modified as. 



This matrix model was studied in [28,29] and similar to the case of // = 0, there is first 
order deconfinement transition at temperature Th- Where Th is given by z{x) — 1. 
There is a discontinuous change in p from zero to a nonzero value. Above Th the 
deconfined phase is preferred and has free energy 0{N'^). Below Th the preferred 
confined phase has free energy 0{1). 

When a small A is turned on, to the quadratic order in A one gets a term of 
the form (Tr(C/)Tr(?7-i)) . One may thus be inclined to propose a matrix model 
corresponding to the black holes with chemical potential as was done earlier with the 
(a, h) model for the black holes zero chemical potential. However, the dependence of 
the coefficients, (a, h) on the chemical potential, though obvious in the A = case, is 
not easy to determine when A 7^ 0. As was mentioned, the lowest correction is O(A^) 
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Figure 9: Plots of the matrix model potential vs. p above jjLc showing the two possible 
situations at finite temperature. At T = the extremal black hole decays into AdS, 
shown as a saddle point at p = 0. 

that involves a three loop computation^^. On the other hand, an alternative approach 
is to write down a model in the strong coupling regime by comparing with gravity. 
Though we do not hope to determine the exact dependence of the parameters on the 
chemical potential, the restrictions on them so that the matrix model reproduces the 
qualitative features of gravity can be inferred. 

Let us begin with the case without a corrections. We will consider the equation 
(31), where now apart from A and T, a and h are also functions of the chemical 
potential 

When comparing with gravity, a(A, T, //) and 6(A, T, p) will be assumed to be 
positive for all values of ji. We have seen in section (2.1) that there exists a critical 
value of the potential (//c) above which there is only one solution. Below this there 
are a maximum of two solutions. The phase structure below the critical potential 
is same as that of the uncharged case. This is reproduced by the (a, h) model when 
a < 1 and 6 > as described in [11]. 

From equation (15), we see that the radius of the small black hole at and above 
the critical potential, becomes zero and negative respectively. The black hole with 
positive radius approaches an extremal limit at T = 0. We know that thermal AdS 
exists for all temperatures and at T = the gauge theory is in the confined phase. We 
thus have two solutions in gravity corresponding to the confined phase. However the 
extremal black hole being nonsupersymmetric will ultimately decay into AdS at zero 
temperature [25,26]. We do not expect to capture the dynamics of this decay with the 

^°To our knowledge this computation is not available in the literature. We hope to return to this 
problem in future. 
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static potential given by the (a, b) model. In fact the zero temperature configuration 
with only the AdS is not continuously connected to the finite temperature (a, b) model 
as we will see below. 

At finite temperature, we always need to introduce an unstable saddle point (a 
maximum) . This follows from the observation that the smooth potential given by the 
(a, b) model which gives two minima (corresponding to thermal AdS and the large 
black hole in gravity) always includes a maximum in between. For /i > /ic this maxi- 
mum corresponding to the unstable black hole ceases to show up in gravity above. In 
the matrix model, we interpret this phenomenon as follows. As n increases beyond 
/^c this unstable saddle point in the matrix model enters the p < 1/2 region from 
p > 1/2 region. Thus the region beyond corresponds to the values of a and b when 
the unstable small black hole of the (a, b) model has already undergone the third 
order Gross- Witten transition [31,32]. This imposes a constraint on the values of a 
and b. Another constraint comes from the fact that the energy of the large black hole 
in this region is lower than thermal AdS for all temperatures. At any finite temper- 
ature, the theory is thus always in the deconfined phase. Both of these constraints 
are satisfied if we set 6 > 2(1 — a). However, in light of these remarks, it is not 
possible to tell whether the unstable saddle point has energy less or more than that 
of AdS or whether it is at p = or away from it. Thus it gives rise to two possible 
scenarios depending on whether the unstable saddle point is at p = or in the region 
0< p < 1/2 as shown in Figure 9. In the following we discuss these two scenarios 
separately . 

(A) Unstable mctximum is at p = {a{p > Pc,T) > 1): In this case we define 
the critical potential Pc by a{pc,T) — 1. Here the unstable saddle point at p = 
does not correspond to thermal AdS but to the unstable configuration not visible in 
gravity as mentioned above. Thermal AdS does not feature in the plot. It has energy 
higher than the black hole. The condition 6 > 2(1 — a) is automatically satisfied as b 
is assumed to be positive. 

(B) Unstable maximum is at < p < 1/2 {a{fi > Pc,T) < 1): In this case the 
saddle point at p = is thermal AdS. The unstable saddle point has energy higher 
than AdS. Since the black hole has energy less than AdS, a and b correspond to val- 
ues above the Hawking- Page transition. Also the saddle point for the maximum is at 
p < 1/2. The critical potential in this case is given by, 6(A, T, pc) = 2 [1 — a(A, T, pc)] 
or by the curve in the parameter space of a and b that gives the Hawking-Page tran- 
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sition, depending on whichever satisfies both the above conditions. 

At this point we recall the plot of [11] showing the various regions in the (a — b) 
parameter space. This is shown in Figure 10. For /i > /ic have the following 
situation: At T = the only solution which is thermal AdS, is given by the region 
below line I. At any finite temperature the parameters jump to values in regions (A) 
or (B). 



IV 




a = 1 



Figure 10: The a — b parameter space for n > The zero temperature region is 
below I and the finite temperature region is beyond III or II. Below I there is only 
one saddle point at p = corresponding to thermal AdS. On II the Hawking Page 
transition occurs. Ill is the Gross- Witten transition line. IV is the line, a = 1 



The matrix model for the fixed charge is obtained by integrating over the chemical 
potential. For this we need to know the exact dependence of a and 6 on in this strong 
coupling regime which we do not have. We note that the possibility (A) incorporates 
a condition only on a. The condition on 6 (6 > 2(1 — a)) is always satisfied in that 
range of a. This implies that we can consistently assume h to be independent of /x. 
However the possibility (B) implies that both a and h should be dependent on /i. 

Equation(33) shows the exact dependence of a on p in the free theory. In this case, 
including only the charged scalars, a(A, T, ji) = [c(A, T) + d{\, T) cos(^)]. If we assume 
this to hold in the strong coupling regime with b independent of fi one can derive a 
model for the canonical (fixed charge) ensemble. This will thus be consistent with 
(A). This fixed charge model was studied in [36]. We will discuss the matrix model 
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corresponding to the Gauss- Bonnet black holes with fixed charge in section (3.5). 
Before this in the next section we study the model with fixed potential including a 
corrections. 



3.4 Including a corrections : Fixed potential 

As mentioned before, once we include d corrections along with nonzero chemical 
potential, the coefficients in the action (32) depend also on a! . The equations of 
motion are of the same form as given in the = case: 

< P < 1, (34) 



4(1 -p) ' 2 

where we have defined 
The potentials that follows from the above action are given by 



= ->l4p' + A3p^-A2p'+(Ai-l/2)p2-llog[2(l-p)] + i , \<P<^. 

As seen from (34) p = is always a solution. The action (32) evaluated at p = 
vanishes. For zero chemical potential and a. this solution corresponds to the thermal 
AdS on the bulk side [11]. 

In the analysis of phase structure we saw as the chemical potential $ and tem- 
perature T vary we arrive different regions having different thermodynamic features. 
Similarly in the matrix model as chemical potential and temperature vary the coeffi- 
cients of the matrix model action vary as well. Thus if we consider a four dimensional 
space corresponding to the four coefficients of (32) it is sectioned into various regions 
which are analogous to the regions in the phase diagram. This is essentially same 
as what we did graphically in Figure 10 for two coefficients. Analogously, here we 
will identify different regions in four dimensional space of the coefficients with various 
ranges of temperature and chemical potential. 

Let us begin with the various regions in the ($^-0;) plane depending on the number 
and nature of the solutions, as discussed in subsection (2.1) (see Fig. 1). In particular, 
there is a line in the ($2-Q!)-plane that separates region with three solutions and that 
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with one solution. Let us consider constraints imposed on the coefficients Ai that 
corresponds to these regions or in other words, regions above and below the critical 
potential /Xc at fixed A. For that purpose it is useful to consider the quadratic poly- 
nomial in p^: /(p2) = (l/p)^F(p). This is given by f{x) = ASA^x^ - 24A3X + 8^2. 
The zeroes of / determine the non-trivial turning points of F. We will see below that 
the parameters at T = are continuously connected to the finite temperature ones 
for p < pc but they are not so for p > pc as we saw in absence of a' correction. 



(A2, A3, A4) 




0.1 0.2 0.3 0.4 0.5 



Figure 11: For < /ic, Potential as function of p for the ranges < p < 1/2 and 1/2 < 
p < 1. The values of {A2, A3, A4 ) used in the plots are given above. Ai = 0.02. 

Sub-critical Potential (p < pc): This corresponds to the region I in Figure 1. In 
this range the behaviour is very similar to the case of zero chemical potential. As is 
evident from the plot in Figure 2(c) obtained on the gravity side we have two charac- 
teristic temperatures, among others. One is T = Tn2 where the stable big black hole 
and unstable intermediate black hole nucleates. Another is T = T/vi > Tjv2 where the 
intermediate unstable black hole combines with stable small black hole and beyond 
that temperature they cease to exist as solutions. This leads to three qualitatively 
different ranges of temperature and in the following we discuss them separately. 

Tn2 < T < Tjvi: Here on the gravity side we have three black holes and so we 
expect at least three solutions of the saddle point equation (34) obtained from the 
Matrix model. The small black hole has a size of the order of a' . In addition, as we 
have already seen in the analysis of p = case [1], we have an (unstable) solution 
which is a maximum of the potential. Since this solution does not appear in the 
gravity that is beyond an analysis of order a' . So we expect these two solutions occur 
generically in the range < p < The other two solutions which are analogue of 
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unstable intermediate black hole and stable big black hole, since they do have bulk 
analogue we expected them to appear generically in the range | < p < 1- Therefore 
we get four solutions for this region. So this region corresponds to the range of the 
coefficients Ai for which (34) has four solutions within the domain of p. Prom (34) we 
see this requires /(p^) = has two positive solutions < p?. < < 1 where p_ is a 
maximum and p+ is a minimum of F{p). That /(p^) = has two solutions requires 

A = 3Al- > 0, > 0, (37) 

where the second condition is to ensure that both the solutions p^ are positive. That 
p± are minimum and maximum of F respectively implies F"[p^) < 0, F"[p+) > 
from which we obtain 

Other conditions that we need to impose so that (34) admits four solutions are 
F{p_) > 1 and F{p^) < 1, which when written in terms of A^s become 



3^4 
A2-A3 



>2A, + mA4^y-^y_ . (39, 

The above conditions (37,38,39, 40) ensure that there exist a minimum (pmin) corre- 
sponding to the small black hole and a maximum {pmax) that does not have a gravity 
analogue. These solutions are in general different from p± and since F is a cubic 
polynomial in p^ the expressions are complicated. In addition, on the gravity side we 
saw that stable small black hole have positive energy. So we expect the solution of 
the saddle point equation that correspond to this stable small black hole should have 
positive energy. So we need at the minimum 

Vl{Pmin)>0 ■ (41) 

whose explicit form is not very illuminating. 

In addition as we said from gravity analysis we expect two solutions should appear, 
generically, in the first half of the domain of p while two other in the second half. 
We consider two halves of the domain of p separately. Let us begin with the range 
< p < |. That we have two solutions of the saddle point equation (34) in this range 
imposes two further conditions: 

P-<1, ^(^)<1- (42) 
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( In this form it is easier to obtain the restriction on the parameters. Otherwise 
we could write < Pmin, Pmax < \ which are more direct but the expressions are 
cumbersome.) These conditions, when written in terms of the parameters, reduce to 
the following two equations: 




(43) 

8 - - ■ 8 - ^''^ 
These conditions (37,38,39,40,41) are real inequalities in four parameters. Each of 
them will give rise to one (or more) codimension 1 wall(s) in the four dimensional 
parameter space described by Ai, A2, A^, A^. Since this involves four parameters it is 
difficult to have a graphical representation of it. As we see each of the inequalities 
corresponds to wall(s) which we cross when we go beyond this range of temperature 
and chemical potential. However, crossing the wall corresponding to (41) will take us 
to some unphysical region as that implies the stable small black hole has energy less 
than that of AdS which implies on the gauge theory side absence of confinement in 
low temperature. The equations (43 and 44) also give rise to codimension one walls 
but a more elaborate analysis is required to understand the correct significance of 
them. 

We expect the other two solutions to appear in the other half namely | < p < 1. 
On this part the situation is pretty much similar to that of (a, 6)-model. This requires 
we have one solution | < pi < 1 such that it satisfies 

where this pi will appear as the analogue of the unstable intermediate black hole 
solution that we got on the gravity side. Actually, this condition is sufficient to ensure 
that there is the analogue of stable big black hole too. This condition (45) along with 
(37,38,39,40) completes the list of necessary and sufficient condition to have four 
solutions of the saddle point equations. As temperature increases, pi will decrease. 
At Gross- Witten temperature [31,32] T = Tg this reaches the lower boundary of this 
region p = |. At T = the restriction on the parameters can be written as 

2^1 + ^ = ^ + 4^ , 3A^ + 8A2>6As. (46) 

o 8 

This corresponds to a third order phase transition. 
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T < Tn2- In this range of temperature the conditions that we obtain on the first 
half of the domain of p remain the same. On the second half of the domain the two 
solutions will merge at T = T^2- This implies in the four parameter space we are on 
the wall that corresponds to (45). The following equations corresponds to Tn2' 

As temperature decreases the value of F{p) will monotonically decreases and there 
will be no solution in the second half. 



T > T/vi: Here we have no solution in the first half of domain of p at T > T/vi- But 
that can happen in two possible ways and we discuss them in the following. One 
possibility is the two extrema of Vi (p) merge at T = Tjvi which implies we cross the 
wall corresponds to (39). In terms of the parameters this means at T — T]\fi 



A2A3 
3Aa 



M. + ie^,[(iiy^-A],.. (48) 



This is a consistent possibility that can occur in the matrix model. Another signif- 
icance of this equation is this corresponds to the bound beyond which we get the 
features of (a, 6)-model and so this region sits in the same universality class of that 
of (a, 6)-model. 

The second possibility is crossing the wall corresponding to (40), where the mini- 
mum of Vi meets the maximum of V2. However, In terms of the coefficients, then, at 
Tjvi we have 



A2A3 



A3 ,2 A,, 



This possibility seems more plausible as it matches with what we saw on the gravity 
side. There at T = T^i the unstable intermediate black hole merges with the stable 
small black hole. At this point we should comment on the relative values of Tm 
and Tg. If Tjvi < Tg this merging occurs before the system reaches it Gross- Witten 
point and therefore there will be no Gross- Witten transition. On the other hand for 
Tjvi > Tg the intermediate solution has already undergone the GW transition and 
becomes a black hole of the order of a'. The other possibility is if T/vi = Tg where the 
small black hole merges with the intermediate black hole exactly at the Gross- Witten 
point or the two walls corresponding to (40) and (43) merge. 
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Figure 12: For ij, > fic, Ai < 0, A < 0: Potential as function of p for the range < /? < 1/2 
and 1/2 < p < 1 .The values of {A2, A3, ) used in the plots are given above. A\ = —0.02. 

Super-critical Potential (/i > fic)'- This region corresponds to the parameter space 
above region I in Figure 1 and therefore should have a single minimum of the po- 
tential. In the matrix model like the fi < fic region here also we have an additional 
unstable maximum of the potential that does not have any bulk analogue and so it 
is beyond the gravity analysis of the a' order. Once again we restrict this in the 
region, p < 1/2 so that this corresponds to some stringy phase as suggested in [11]. 
Therefore, similar to the analysis for a' = ( and unlike the p < fic case ) we consider 
two possible scenarios depending on whether the unstable maximum lies at p = or 
it lies away from p = 0. In the following we discuss the two scenarios separately. 

(A) Unstable maximum is at p = {Ai < 0): If we have [Ai < 0) the unstable 
maximum will always be at p = 0. In addition, in the matrix model we expect to 
have one and only one stable minimum for p > (The energy of this black hole could 
be greater or less than that of thermal AdS depending on the values of p and A (see 
Figure 3). The condition Ai < alone does not ensure this and we need to impose 
additional constraints on A2, A^ and A4. We note that there are two ways in which 
one can ensure that there is only one stable black hole solution. We discuss them 
separately in the following: 

One possibility is, at p = pc the three solutions merge into one at some value of 
p = P4. > 0. A similar merging occurs on the gravity side when the chemical potential 
reaches its critical value. Since the corresponding solution is manifested in the gravity 
limit we expect that p+ > 1/2. Then, at p = pc, apart from the condition Ai = 0, 
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A2, A3 and A4 satisfy, 

V'{p+) = ; V"{p+) = ; V'" {p+) = for ^ < p+ < 1. (50) 

Thus we obtain the parametric solutions of A2{p+), A^i^p^) and A^i^p^) from (50) 
which corresponds to the hne separating the number of solutions in gravity in Figure 
1. We do not write these parametric equations here as they are not compact enough. 
For, p+ = 0.8 we have, Ai = 0, A2 = 0.583, A3 = 1.481, A4 = 1.293. 

Examining the equation of motion (34), one can find an indirect way to obtain 
a (slightly more general) constraint in compact form. This ensures there is exactly 
one solution at some point p > in addition to the one at p = 0. However, that the 
three solutions merge at this point is not guaranteed. If we relax the first inequality 
of (37) by making the discriminant A negative so that there is no turning point of F 
then we are left with only a single solution. This can be thought of as crossing the 
wall corresponds to the first condition of (37) with p — pc- The merging of the three 
solutions at this point appears as a special case of it and so that condition is more 
ramified. We have plotted the associated potentials in Figure 12. 

The other possibility that one can consider to obtain a single minimum is setting 
A^Ip — pc) — and A4{p > pc) < 0. This amounts to relaxing the second inequality 
of (37) so that one of the turning point becomes imaginary. This corresponds to iden- 
tifying p = Pc with the second condition of (37). For this choice, as p increases beyond 
its critical value, one of the turning point comes down to p = and then disappears. 
This possibility is not continuously connected with the above mentioned constraint, 
which gives rise to merging of three solutions. However, this possibility can give a sim- 
ple form for the fixed charge case. We have plotted the related potentials in Figure 13. 

(B) Unstable maximum lies at0<p<l/2 {Ai > 0): In this scenario for p > pc 
the unstable maximum is in the range < p < 1/2. The usual saddle point p = is 
also there but this time it corresponds to the thermal AdS. In addition, there exists 
only one minimum. (The energy of this minimum is either greater or less than that 
of AdS at low temperatures. If the energy is greater than that of AdS energy, at high 
temperature, the black hole undergoes Hawking Page transition.) In order to make 
sure that this is the one and the only one minimum, once again we identify pc to be 
the limit where three saddle points merge. However unlike the previous case where 
the small maximum was at p = 0, here it lies in0<p<l/2. This configuration 
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Figure 13: For jx > ^c, < 0, ^4 < 0: Potential as function of p for the range < p < 1/2. 
The values of {A2, A^, A4 ) used in the plots are given above. Ai = —0.1. 

satisfies, 



V'{p+)=0 ; V"{p+) = ; V'" ip+) = for ^ < P+ < I- 

and V'{p-) = for < p_ < ^. 



(51) 




Figure 14: For p > pc and Ai > 0: Potential as function of p for the range < /? < 1/2 
and 1/2 < p < 1. The values of {A2, A3, A4 ) used in the plots are given above. Ai = 0.298. 

Apart from being functions of p+ the Aj's are now also functions of p_. With 
p_ = 0.5 and p+ = 0.8 we get, Ai = 0.418, A2 = 1.561, A3 = 2.501 and A4 = 1.691. 
In this case the energy of the merging point (p+) is positive. This gives rise to a stable 
minimum with positive energy as we move beyond p^c- We now vary the parameters 
with increasing temperature and see that the minimum crosses zero corresponding to 
the Hawking-Page transition in the bulk. The plots for p, > p,c are shown in Figure 
14. 
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One can similarly derive the parameters for which the energy of the merging 
point is less than zero. This is given by, p_ = 0.4 and p+ = 0.72 which gives 
Ai = 0.006, A2 = 0.06, A3 = 0.291 and = 0.535. Note that for p_ = we get 
{Ai = 0), which is the possibility (A). As we move away from the critical potential, 
this gives rise to a minimum with energy less than AdS. With further variation of the 
parameters this saddle point goes deeper as shown in Figure 15. This can be mapped 
to the variation of the minimum with temperature as it happens for the black hole 
in gravity. 
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Figure 15: For > (pcAi > 0: Potential as function of p for the range < p < 1/2 and 
1/2 < /O < 1. The values of {A2, A3, A^ ) used in the plots are given above. A\ = 0.002. 

In the above discussion we consider, on the matrix model side, the various possibil- 
ities, which can reproduce the bulk behaviour that we obtained in the analysis on the 
gravity side. The different possibilities gives rise to different ranges of the parameters, 
of which some are mutually exclusive. With this amount of input it is not possible 
to decide which one is the correct behaviour. We need the explicit dependence of 
Aj's on p in the strong coupling regime. Perhaps a weak coupling calculation of the 
associated terms in the model can serve as a clue. 

3.5 Including a corrections : Fixed charge 

We have seen in the last section that when restricted to various values of the param- 
eters Ai, the matrix model (32) reproduces the features of gravity below and above 
the critical potential, In general thus Aj's will depend on the chemical potential. 
In the study of the matrix model for fixed charge, the explicit knowledge of the de- 
pendence of Aj's are necessary. Out of the cases studied in the earlier section, there is 
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one that allows us to construct a toy model for the fixed charge consistently. This is 
given by the second possibility of (A), when both Ai and A4 become negative above 
the critical potential. 

In the following, we will assume the result from the zero coupling to be valid in this 
strong couphng regime. In the free theory, Ai is given by, (1 — 2Ai) — [c + dcos(//)] 
when matter fields consist of only charged scalars. Since we are only interested in the 
qualitative features of this model, for simplicity we will take, A4 — a4(A, T) cos(//) 
with A2 and As functions only of A and T. The action with this dependence on the 
potential is. 



ai cos(/x)p^ - Asp^ + Aap^ + 



' c + dcos{p) 



P 



(52) 



The partition function for the canonical ensemble is obtained by integrating over 
the chemical potential. Hence following [36], we write. 



Z{Q,T,X) = j diJ,exp{-iij,Q) J dpexp{-N'^V{p, p)}, 
where, the potential is. 



(53) 



Vip) = -^^e//(p^) + 



with. 



■^5e,,(p^)-llog[2(l-p)] + - 



< p < 
1 



2' 
1 



<P<1, 



(54) 



Seffip') = N^[cp^ + 2A2p'-2Asp']+log[lQ(N\2aip' + dp^))]. (55) 

Here Iq{x) is a Bessel Function. We are interested in the the large N hmit. This is 
obtained by keeping in mind that Q"^ ~ C(iV^) so that Q'^ = N'^q where q ~ 0{1). 
The resulting effective action in the large N hmit is thus. 



+ N^q 



cp" + 2A2p^ - 2Asp' 

' p^{d + 2aip^] 



+ 



1 

2\ 2 



log 



1 



1 + {1 + 



p4(c;+2a4p6)2 



(56) 



32 



V(P) 




Figure 16: Coloured curves: Potential for various temperatures corresponding to the 
variations of Ai, A2, A3, a, c and for q < Qc Black curve: Potential for q > q^. 

Defining, F{p) = dS{p^)/dp^ , the equation of motion is given by (34). 

The phase structure is qualitatively the same as when there are no a corrections. 
Unlike the fixed potential thermal AdS is not a solution. It is easily seen from (56) 
that p = is not a solution of the equation of motion. We have seen in section 2.2, 
that there exists a critical charge for a fixed a beyond which there is only one solution. 
Below this critical charge there are a maximum of three solutions. See equation (23) 
and Figure 5. In the matrix model, this critical charge (gc) is given by the merging 
of the three saddle points of (54). 

Figure 16 shows the thermal history for the canonical ensemble. The coloured 
curves are for q < qc for some fixed A. This corresponds to region / in Figure 5. At 
low temperatures, there is a single minimum in the range < p < 1/2 (curves in 
pink). This corresponds to the small black hole. Two new saddle points appear at 
T = T3 with the new minimum in the range 1/2 < p < 1 (shown in blue) . These 
correspond to the intermediate (unstable) and large black holes. The small black 
hole has lower energy upto Tc (the green curve), beyond which the large black hole 
phase is favoured (yellow). Thus there is a phase transition at corresponding to the 
one in the bulk. The small and the intermediate black holes disappear at Ti and at 
high temperature, the large black hole (the phase with the saddle point in the range 
1/2 < p < 1) remains (shown in purple). 

The black curve in Figure 16 shows the single saddle point when q > q^. This 
corresponds to the single large black hole that exists above the critical charge for 
fixed a, in the region II of Figure 5. 
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